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Abstract 

We study Topological Defects (TD) in hidden (mirror) matter as possible 
sources of ultra-high energy neutrinos. The hidden/mirror and ordinary mat- 
ter are assumed to interact very weakly through gravity or superheavy par- 
ticles. An inflationary scenario is outlined in which superheavy defects are 
formed in hidden/mirror matter (and not in ordinary matter), and at the 
same time the density of mirror matter produced at the end of inflation is 
much smaller than that of ordinary matter. Superheavy particles produced by 
hidden-sector TD and the products of their decays are all sterile in our world. 
Only mirror neutrinos oscillate into ordinary neutrinos. We show that oscil- 
lations with maximal mixing of neutrinos from both worlds are possible and 
that values of Am 2 , needed for for solution of solar-neutrino and atmospheric- 
neutrino problems, allow the oscillation of ultra- high energy neutrinos on a 
timescale of the age of the Universe. A model of mass-degenerate visible and 
mirror neutrinos with maximal mixing is constructed. Constraints on UHE 
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neutrino fluxes are obtained. The estimated fluxes can be 3 orders of mag- 
nitude higher than those from ordinary matter. Detection of these fluxes is 
briefly discussed. 

I. INTRODUCTION 

A hidden sector of mirror particles was first suggested by Lee and Yang [|IJ in 1956 to save 
the conservation of parity in the whole enlarged particle space. This concept was further 
discussed and developed in Ref . [Q] . Later the idea of two weakly interacting sectors, visible 
and hidden, found interesting phenomenological applications and development [0,0. It has 
been boosted in 1980s by superstring theories with E s x E' s symmetry. The particle content 
and symmetry of interactions in each of the E$ groups are identical, and thus the mirror 
world has naturally emerged. 

The most recent reincarnation of hidden-sector models is in the context of D-branes 
In this approach, light particles are associated with the endpoints of open strings which are 
attached to D-branes. Ordinary and hidden-sector particles live on different branes which 
are embedded in a higher-dimensional compactified space. 

How do the ordinary and hidden sectors communicate with each other? 
Most naturally they interact gravitationaly. This possibility is employed in J7|. More gen- 
erally, and this is also discussed in |7[] and ||, ordinary and mirror matter can also interact 
through the exchange of superheavy gauge particles. In the D-brane context, in some models 
the interaction between different branes occurs only through the exchange of closed strings 
(gravitons), while other models (in which the two branes are embedded in a brane of higher 
dimensionality) allow for a gauge boson mediated interaction. 

In the case of gravitational coupling of the two worlds, one should expect dimension 5 
gravitational interaction scaled by the Planck mass m p i; in the case of superheavy gauge 
bosons this scale might be A < m v [. 
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where ip and H are respectively the lepton and Higgs SU(2) doublets, with mirror fields 
denoted by primes. Eq. (|1|) provides mixing of ordinary and mirror neutrinos and neutrino 
masses ||, |I0 |. Because of the smallness of the neutrino masses, this is the most visible 



physical effect caused by the gravitational interaction (A = m p i) of particles from the two 
worlds. 

In principle, there could be other ways of communication. For example, (see |12|]) one 
can add to the Lagrangian a Higgs potential term A000'0' an d a gauge boson kinetic mixing 
term hF^F^, with A and h being new coupling constants. These terms have potential 
problems. 

The discrete P-symmetry that interchanges the two worlds can be spontaneously broken. 
In this case, the coupling constants, the Higgs potentials and expectation values, and even 
the symmetry breaking patterns will generally differ from one world to the other. The 
breaking of P-symmetry can be implemented by giving a non-zero vev to a spin-0 field 
which transforms as a singlet under the gauge groups in both sectors ans as a pseudoscalar 
under the P-transformation: P0 = —0. Models of this type have been studied in Refs. 
|7|jn|| ; we shall refer to them as asymmetric hidden sector models.^ 



A model with an unbroken P-symmetry has been developed in Ref. [|1J|. The P- 
transformation in this model turns the left-handed ordinary fermions into right-handed mir- 
ror fermions. The masses and couplings of ordinary and mirror particles are identical, and 
hence the term "mirror matter" is more justified in this case. The EW Higgs fields in both 
sectors are also parity partners and have equal vev's. Mixing of neutrinos from different 



1 One can, of course, consider models without P-symmetry, in which the ordinary and mirror 
sectors are not symmetric to begin with. For our purposes, such models are essentially the same 
as models with a spontaneously broken P-symmetry. 



worlds (taken as an ad hoc terms in the Lagrangian) provides a channel of communication 
between the worlds: ordinary neutrinos can oscillate into mirror neutrinos, which play the 



role of sterile neutrinos. We shall refer to models of this type (with model [12] as the most 
elaborated example) as symmetric hidden sector model. 

Although the term "mirror sector" suggests that it is related to the ordinary sector by a 
reflection transformation, most of the discussion in this paper is applicable to a more general 
class of hidden-sector models. In what follows, we shall use the terms "mirror sector" and 
"hidden sector" interchangeably. 

Mirror neutrinos are the most natural candidates for sterile neutrinos often considered 
now for explaining the oscillations of solar and atmospheric neutrinos. In this respect mirror 
neutrinos have been studied in Refs. |7|JT2|. 

Mirror matter has cosmological consequences which result in model restrictions. 

In symmetric hidden sector models, the number of massless and light particles is doubled 
in comparison with ordinary matter, and this case is excluded by cosmological nucleosynthe- 
sis, if the temperatures of mirror and ordinary matter are the same. (More generally, in any 
mirror model the effective number of light degrees of freedom is larger than for the ordinary 
matter, and this number is restricted by nucleosynthesis.) 

One way to suppress the light degrees of freedom is to diminish the temperature of the 
mirror matter in the Universe [[12| , |ll~| (see further discussion in Section II). This reduces the 
number density of mirror photons in a sraightforward way, while the situation with sterile 
neutrinos is more delicate [O]. Even if the initial density of mirror neutrinos is negligible, 
they reappear again and may be brought to equilibrium at nucleosynthesis epoch due to 
oscillations between ordinary and mirror neutrinos |l4j]. Nucleosynthesis constraint bounds 
the allowed neutrino properties in the parameter space (Am 2 , sin 2 29), where 9 is the mixing 



angle jlSj . It is clear that the larger the mixing angle is, the smaller are the allowed values 
of Am 2 . Electron neutrinos impose the strongest limit on Am 2 , because v' <-> u e oscillation 
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influences nucleosynthesis not only through the rate of the cosmological expansion, but also 
due to the rate of n <-> p conversion. 

A cruicial assumption involved in deriving the bounds described above is that the relic 
lepton asymmetry is absent. In the presense of a large lepton asymmetry, L > 5 • 10~ 5 , 
the potential for active neutrinos results in a small mixing angle between active and sterile 
(mirror) neutrinos, and thus in a weak oscillation between these components at temperatures 
relevant for nucleosynthesis [13], [16], [17]. Therefore, in this case v' v oscillations and 
induced by them the additional number of degrees of freedom are suppressed. 

A cosmological origin for the temperature difference between mirror and ordinary matter 
in the universe has been already considered in Refs. [|J , jjnj . The asymmetry is generated as 
a result of different rates of inflaton decay to ordinary and mirror matter due to P-symmetry 
breaking. 

In Ref. jnj mirror neutrinos and baryons were considered as dark matter particles. The 
problem of structure formation with mirror matter and other astrophysical implications were 



studied in Refs. 11 , [IS 1, [PL 



In this paper we shall study mirror matter in the Universe as a source of ultra-high 
energy neutrinos. As concrete sources we shall consider mirror Topological Defects. They 
produce high energy mirror neutrinos in the usual way: through production and decay 
of superheavy mirror X-particles. Then high energy mirror neutrinos oscillate into ordinary 
neutrinos, while the other products of decay of mirror X-particles remain in the mirror world, 
being invisible in the ordinary matter. These sources give an ideal example of "hidden 



2 In the scenario studied in this paper only muon neutrinos from ordinary and mirror worlds are 
maximally mixed, while the other two types of neutrinos might have very small mixing. Therefore, 
even in case L = 0, this scenario gives only one extra neutrino and thus marginally survives the 
nucleosynthesis restrictions 
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neutrino sources" |20|. High energy neutrino radiation from ordinary sources is usually 

accompanied by other radiations, most notably by high energy gamma rays. Even in cases 

when high energy photons are absorbed in the source, their energy is partly transformed into 

low energy photon radiation: X-rays or thermal radiation. The fluxes of these radiations 

impose an upper bound on the high energy neutrino flux. For sources transparent for HE 

gamma radiation, in particular for Topological Defects (TD), the upper limit on diffuse 

neutrino flux is imposed by the cascade e-m radiation PHI , IPOfl . In all cases (e.g. decay 

of X-particles, pp and pj interaction) neutrinos and photons are produced at the decays 

of pions. Colliding with microwave photons, high energy photons and electrons produce 

e-m cascade with most of the energy being in the observed 100 MeV — 10 GeV energy 

range. The energy density of this cascade radiation should not exceed, according to EGRET 

observations, u) cas ~ 1 — 2 • 10~ 6 eV/cm 3 . Introducing the neutrino energy density for 

neutrinos with individual energies higher than E, u> u (> E), it is easy to derive the following 

chain of inequalities (from left to right): 

Att r 00 Att r°° Att 

u cas > lu u (> E) = — EI v {E)dE > —E / I v (E)dE = —EI U (> E), (2) 

C JE C JE C 

An upper bound on the integral diffuse neutrino flux immediately follows From Eq.(H): 

u> E) < s-r- (3) 

The upper bound given by Eq.(|3|) is the most restrictive one for diffuse neutrino fluxes 
produced by "ordinary" TD, AGN, gamma-ray bursts, etc. The neutrino flux from hidden- 
sector TD is free from this bound, though in Section V we shall present another (weaker) 
upper limit on high-energy neutrino flux from hidden-sector TD. 

Another problem where UHE neutrinos from mirror TD can be helpful is the production 
of Ultra High Energy Cosmic Rays (UHECR) by resonant neutrinos. 



The signature of extragalactic UHECR, the Greisen-Zatsepin-Kuzmin (GZK) cutoff |22 



is not found in observations [23]. Three ideas have been suggested to explain the absence of 



the cutoff: (i) signal carriers are not absorbed on the microwave radiation, (ii) the sources 
form a compact group around our Galaxy with a linear size smaller than the GZK absorption 
lenght, and (iii) the sources are distributed uniformly in the Universe, but the target particles 
on which UHECR are produced by the signal carriers form a compact object. 

The case (iii) is realized [^4 with the help of UHE neutrinos (signal carriers) interacting 
with relic neutrinos which have an enhanced density in some Large Scale Structure near us. 
One of the problems with this idea is that it is not clear how one can arrange a large flux 
of UHE neutrinos. The origin of this flux is left unspecified in all publications (e.g. P5f| , 
26]), [27] that we are aware of. The very high resonant neutrino energy, Eq ~ 4 ■ 10 12 GeV 



for m v ~ 1 eV, implies a top-down scenario, but the fluxes of neutrinos in such scenarios 
are limited by the cascade constraints and are too small to produce the observed flux of 
UHECR. Mirror TDs, which evade the cascade constraints, can in principle produce the 
desirable UHE neutrino flux. This problem is addressed in Section VI. 
The ouline of this paper is as follows. 

In Section II we review inflationary scenarios for the mirror Universe, elaborating in 
particular a two-inflaton scenario for symmetric HS models. Formation of hidden-sector 
TD is discussed in Section III. Fluxes of High Energy Neutrinos from TD are calculated 
in Section IV for the case of necklaces as an example. An upper bound on the neutrino 
flux from hidden-sector TD is obtained in Section V. UHECR from resonant neutrinos are 
studied in Section VI. The main results of this paper are summarized in Section VII. In 
Appendix our model for mass-degenerate neutrinos with maximal mixing is described. 

Throughout the paper we shall use following abbreviations: HE and UHE for High Energy 
and Ultra High Energy, respectivey, UHECR - for Ultra High Energy Cosmic Rays, TD - 
for Topological Defects, HS - for Hidden Sector (including mirror sector), LSS - for a Large 
Scale Structure in the universe, LG - for the Local Group of galaxies, LS -for the Local 
Supercluster of galaxies. 
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II. INFLATIONARY SCENARIOS FOR HIDDEN-SECTOR UNIVERSE 



As was discussed in the Introduction, there are several dangers to be watched for in 
models with mirror matter. The main one is a possible conflict with the standard nucle- 
osynthesis. If the mirror sector contains a massless photon and three light neutrinos, and 
the temperatures of the two worlds are the same, then the density of mirror matter at the 
time of nucleosynthesis is unacceptably high (it amounts to five extra neutrino species). 

Two conditions are necessary to overcome this problem (see Introduction): The mirror 
matter must have a lower temperature T" ^ 0.5T and lepton asymmetry is needed to suppress 
excessive production of mirror neutrinos through oscillation of ordinary neutrinos. Here, we 
shall discuss some inflationary scenarios that can naturally lead to a temperature difference 
between the two worlds. 

The condition T' % 0.5T can be easily satisfied in asymmetric hidden sector models. For 
example, we could have a single inflaton field which transforms as a scalar (or pseudoscalar 
0) under the P-transformation. Since the symmetry between the two sectors is broken, the 
field will generally have different couplings to particles in different sectors. It will then 
decay into ordinary and mirror particles at different rates, and the reheating temperature in 
the mirror matter can be lower HJIT| . 

This scenario would not work in symmetric hidden sector models: the inflaton would 
then have identical couplings to both sectors, and the two reheating temperatures would be 
the same. This problem can be addressed in the following two-inflaton scenario (see also 

!)■ 

Let us consider two inflaton fields, and <f)', with belonging to the visible sector and 
0' to the mirror sector. During inflation, both inflatons roll down towards the minima of 
their respective potentials. Inflation ends when both of them have reached their minima. 
An important point is that the evolution of and 0' need not be synchronized. The inflaton 
dynamics is influenced by quantum fluctuations which cause inflation to end at different 
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times in different regions of space p8| . In our case there are two inffatons, their fluctuations 
are uncorrelated, and one expects them to reach their minima at different times, even in 
the same spatial regions. In regions where 0' reaches minimum first, any mirror particles 
produced due to its oscillations are diluted by the remaining inflation driven by the field 
0. By the time when the energy of thermalizes, the density of mirror matter will then 
be very small, so that T ^> T". Note that the (co-moving) coherence length of the inflaton 
fields should be much greater than the present horizon, so we expect very large (super- 
horizon) regions of the universe to be dominated by ordinary matter, and similarly large 
regions dominated by mirror matter, with relatively tiny boundary regions where both kinds 
of matter are present in comparable amounts. It is very unlikely for us to find ourselves in 
one of such rare regions. 

A quantitative analysis shows that this de-synchronization picture may or may not apply, 
depending on the form of the inflaton potential £7(0). We shall see, however, that it does 
apply for the simplest choice of the potential, 

C/(0) = my/2. (4) 

In "chaotic" inflation scenario, the inflatons roll from very large values of towards 
= 0. The initial values of and 0' are large and uncorrelated. At very large 0, > <p q , the 
dynamics of are dominated by quantum fluctuations. The boundary g of this "quantum 
diffusion" regime is determined by condition dJJ/dcj) ~ H 3 . At 0* < < <p q , where 0* ~ 
rripi, (and 0') evolve in a slow-roll regime described by the equations 

3^=-^, (5) 



dUU') 

3H<P> = -%±, (6) 



H 2 = ^[U{<P) + U{<P% (7) 
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Will the "incidental" initial ratio 4>'/<p be conserved during the slow-roll evolution? For 
the potential @) the answer is "yes" . Indeed, the integration of Eqs. © - (0) results in 

0'/0 = const. (8) 

Suppose that 0' at the time when begins its slow roll (0 ~ <p q ~ m pl ), so 

that we can neglect U{<p') in (|7]). The characteristic value of 0' is then determined by quantum 
fluctuations about 0' = and can be found from [^9] U(4>') ~ H^, where H q ~ (m^m^) 1 / 2 is 
the expansion rate at ~ 9 . This gives 0' ~ m p ; and 0'/0 ~ (m^,/m p ;) 1//2 . This small ratio 
of 4>'/4> is preserved all the way to the end of inflation. In this type of models, the universe is 
divided into super-horizon regions dominated by ordinary matter and equally large regions 
dominated by mirror matter. 

For a different choice of the potential, e.g. 

f/(0) = A^(0/m pi ) n , (9) 

with n > 2, the fields and 0' do get synchronized at late stages of the evolution. (For 
n = 4, this was shown in Ref. |30[| .) In this case, integration of the slow-roll equations (|5|) 

-(ED g ives 

0- n+2 - 0'~ n+2 = const. (10) 

At the onset of the slow roll of 0, when ~ <p q ~ A _1 / 4 m p ;, the typical value of 0' is 
0' ~ A^ 2 n ^ n ^ n+2 ^m p i. With < 1, this is much smaller than 9 but still much larger 
than 0*. Now, consider the solution (|T0| ) with these initial values of and 0'. By the end 
of inflation, both fields get much smaller than their initial values, so that the constant on 
the right-hand side of (|T0|) becomes unimportant, and we have 0' ~ 0. Thus, and 0' get 
synchronized by the end of inflation, even if they were not initially. We conclude that models 
with a power-law potential @ and n > 2 give equal densities of mirror and ordinary matter. 
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Coming back to our basic ordinary-matter dominated scenario, we can give another 
example of a model with segregated mirror and ordinary matter is a two-inflaton model 
where inflation occurs at a metastable minimum of the inflaton potential. The highest rate 
of inflation is achieved in the false vacuum state where both and are at the minima of 
their respective potentials. This state decays through nucleation of two types of bubbles. In 
bubbles of the first type, the field <ft' tunnels through the barrier and starts rolling down its 
potential, while remains in the false vacuum. As 0' rolls to the bottom of the potential and 
decays into mirror particles, inflation continues in the interior of the bubble. Mirror particles 
are quickly diluted away, and the bubble interior is filled with inflating false vacuum of the 
"ordinary" inflaton 0. This vacuum will in turn decay through nucleation of bubbles of the 
second type (with tunneling and 0' remaining unchanged). If the potential is sufficiently 
flat, the roll down of is accompanied by additional inflation, and after the 0-field decay, 
the interiors of these secondary bubbles will become dominated by ordinary matter. Quite 
similarly, nucleation of (//-bubbles inside 0-bubbles results in regions dominated by mirror 
matter. 



III. HIDDEN-SECTOR TOPOLOGICAL DEFECTS 

Apart from the general considerations which apply to any model with a hidden sector 
(HS), we have to address some additional issues specific to HS defects. First we have to 
arrange for these defects to form. And second, we have to avoid the formation of similar 
defects in our sector, since otherwise we would get unacceptably large fluxes of ordinary 
UHE particles, and the cascade bound would be violated. 

Once again, these conditions can be easily satisfied in asymmetric HS models, where the 
symmetry breaking scales and even the symmetry breaking patterns may be different in the 



two worlds. Perhaps the simplest possibility is the model of Ref.( llj) with one inflaton 



and asymmetric reheating. HS defects can be formed in a usual symmetry-breaking phase 
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transition after inflation. We only have to arrange for the corresponding phase transition in 
the ordinary matter to occur at a lower energy scale or not to occur at all. 

In symmetric HS models, the two sectors have identical physics and identical defect 
solutions. It does not follow, however, that ordinary and HS defects should be present in 
the universe in equal numbers. The density of defects is determined by the cosmological 
evolution, which can be different for the two sectors. For example, in the two-inflaton model 
discussed above, mirror matter is completely inflated away, and if defects were formed in 
phase transitions after inflation, we would expect to have ordinary defects but no HS defects. 
For our purposes, however, we need the opposite situation: HS defects and no ordinary 
defects. 

This can be arranged if defects are formed in a curvature- driven phase transition during 
inflation. 

As an illustration we shall consider a toy model of a spontaneously broken 577(2) sym- 
metry. We introduce a Higgs triplet x — (Xi? X2, X3) with a potential 

v(x) = \h\x\ 2 -v 2 ) 2 - (ii) 

When x acquires an expectation value, SU(2) is broken to U(l) and monopoles are formed. 
Suppose now that x is coupled to the inflaton (ft, 

VM = -\g<t> 2 \x\\ (12) 

and has a non-minimal coupling to spacetime curvature, 

3 In models discussed in Refs. PJTl[|, the P-symmetry breaking is assumed to be at the electroweak 
scale, rjp ~ t\ew ~ 10 2 GeF, and the two worlds have nearly identical physics above t]ew- But this 
needs not be the case. The physics of mirror and ordinary defects will be different in models where 
r]p is greater than the symmetry breaking scale of the defects. 
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Vr( X ) = ^Ixl 2 - (13) 

The mirror field x' nas identical couplings to the mirror inflaton 0' and to the curvature. 

We shall assume "chaotic" inflation with a quadratic inflaton potential (f|). After the 
mirror inflaton cf)' reached its minimum, the effective mass of the mirror field \' is 

m 2 x , = -At/ 2 + £R, (14) 

and the curvature is given by 

R « 12H 2 « 167r(m /m p ;) 2 2 - (15) 
We see that above the critical curvature 

Rx> = w/t (is) 

the expectation value of x' vanishes and the SU(2)' symmetry is restored. Thus, even in the 
absence of mirror matter, as the curvature gradually decreases during inflation, we can have 
symmetry-breaking phase transitions accompanied by the formation of HS defects. If these 
curvature-driven phase transitions occur sufficiently close to the end of inflation, the defects 
will not be completely inflated away and can serve as sources of UHE neutrinos. 
The effective mass of the field x i n t ne ordinary sector is 

m\ = -Ar/ 2 - gtf + £R = -At? 2 - g[l - 16tt(£/ ' g^/m^ 2 . (17) 

We see that m 2 < for all <fi, provided that 

A 

lQ^lg){m^/m pl ) 2 < 1. (18) 

In this case no ordinary defects are formed during the whole slow roll period of (and any 
defects formed prior to that are completely diluted away) 

One final condition that has to be checked is that the temperature at reheating is not 
so high that the symmetry gets restored, since otherwise "ordinary" defects will be formed 
again in a subsequent phase transition. All the above conditions can be satisfied without 
fine-tuning. 
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IV. FLUXES OF HIGH-ENERGY NEUTRINOS FROM HS DEFECTS 



We shall consider necklaces as a specific example of TD. Necklaces are hybrid TDs 
formed by monopoles (M) and antimonopoles (M), each being attached to two strings. 
The monopole mass m and the mass per unit length of string /x are determined by the 
corresponding symmetry breaking scales , rj s and r) m , 

m ~ ATrr) m /e, fJ, ~ (19) 

where e is the gauge coupling. The evolution of necklaces depends on the parameter 

r = mj fid (20) 

which gives the ratio of the monopole mass to the average mass of string between two 
monopoles (d is the average string length between monopoles). It cannot exceed r max ~ 
T] m /r) s . As it is argued in ref. P^] , necklaces might evolve towards a scaling solution with 
a constant r ^> 1, possibly approaching r ~ r max . Monopoles and antimonopoles trapped 
in the necklaces inevitably annihilate in the end, producing superheavy Higgs and gauge 
bosons (X particles) of mass mx ~ er] m . The rate of X-particle production per unit volume 
and time is 

ii x ~ r 2 fi/t 3 m x (21) 
From the relations above it is easy to see that 

C-^-§r— V 2 S10. (22, 

X 



ni 



High energy neutrinos are produced in the chain of X-decays via pions. For simplicity 
we assume that pions (of all charges) are produced with a power-law spectrum 

D v (x, m x ) = 4(2 - p)2~ p x- p (23) 
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where x = E n /mx is a fraction of energy taken away by a pion, and for p we shall use a 
value between 1.3 and 1.5 , which bound the realistic QCD spectrum of pions. The spectrum 

1 /2 

(j23|) is normalized so that J xD n (x,mx)dx = 1. 

The number of neutrinos with energy E v from the decay of one X-particle is given by 

4 r 1 / 2 dx 

N{E U ) = / —D n {x,m x ). (24) 

ITlx J2E v /m x X 

The diffuse flux of mirror neutrinos v[ (where i = e and fi , antineutrinos are not included) 
is 

c nxjto) [™ f 1 / 2 dx 
Iv i (E) = — / dz — D 7T (x,m x ) (25) 

HmxtLQJa J2Ev{l+z)/m x x 

Finally, we obtain the diffuse flux of ordinary neutrinos Vi taking into account v[ — > z/j 
oscillation with averaged probability P osc ~ 1/2: 

ue) = ±^(^\-'(1^ Pck (26) 

Anm x H \mx/2j p(p - 1) 
This expression can be written in more compact form: 

c 1 (E v 



with ( given by Eq.(|22|) and 



kpC-r* — . ( 27 ) 



_Q(2-p)2~ 2p 

p(p- 1) 

For P osc = 1/2 and p = 1.5 fc p = 1/4 and A; p = 0.89 for p = 1.3. 

The neutrino flux in Eq.(^) can be very large. For example, with r 2 /i ~ 0.1m x , nix ~ 
10 16 GeV, p = 1.5 and P osc = 1/2, one obtains at E ~ 10 11 Gel/, E 3 ^ ~ 1 ■ 
10 28 eV 2 m~ 2 s~ x sr~ x , i.e. a flux three orders of magnitude larger than predicted from ordi- 
nary sources under most optimistic assumptions. Since we have used a very rough power-law 
approximation for the spectrum, it may be better to illustrate the enhancement of the flux 
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by comparing Eq. ( ^6|) to the maximum allowed flux from ordinary sources with the same 

power-law spectrum. 

Suppose that some unidentified ordinary-matter sources produce superheavy X-particles, 
which decay producing high energy pions. Then the energy density of neutrinos uo v and the 
cascade density uj cas are equal. Using this fact we immediately obtain an upper limit on the 
neutrino flux as 



The ratio of the two fluxes [Eq.(^) and Eq. ([29|) l is given basically by the value m 2 x /(uj cas tl) 
and it is ~ 2 • 10 4 for nix ~ 1 • 10 16 GeV and for observationally allowed cascade density 
u cas ~ 2 ■ 10~ 6 eV/cm 3 . 



All particles from the QCD cascades produced by decays of mirror X-particles are sterile 
in our world. Only mirror neutrinos can oscillate into ordinary ones. An upper bound on the 
neutrino flux is given by the resonant interaction of UHE neutrinos with relic cosmological 
neutrinos, v + v — > Z° — > pions. Electrons and photons from the decay of pions initiate e-m 
cascades on the microwave radiation. Reactions v + v — > Z° — > / + I, where I = e, /i, r, also 
contribute to the cascade. The calculated cascade energy density u cas must be smaller than 
the energy density u° bs observed (e.g. by EGRET) in the extragalactic diffuse radiation. 

Let us first derive a convenient formula for the rate of resonant events. 

The resonant neutrino energy E and the resonant cross-section u(E) for u + u — ► Z° — > / 
(/ is an arbitrary final state) are given by 




(29) 



V. E-M CASCADE RESTRICTIONS 




(30) 



°vj( E c) 



12tt Y u Y f 
ri^(E c -m z y + Tl/A 



(31) 



16 



Here, mz is the mass of Z°-boson, E c is the center-of-mass energy, T u , Tf and T t are the 
widths of Z° decay to neutrinos, to an arbitrary final state / and the total width, respectively. 



In Eq.(31) we took into account that only one chiral component of neutrino takes part in 
the interaction. [] 

The rate of Z° production per unit volume due to collisions of high energy flavor neutrino 
Vi with target antineutrino i/j is given by 

h z = J I Ui {E)a{E)dE = 4ira t n p £I Ui {E )E , (32) 

where 

a t = ASixf v G F = 1.29 ■ 1(T 32 cm 2 , (33) 

Gp is the Fermi constant, and here and below f s = T s /GFm 3 z , with T s being the width of 
the channel. Numerically, f v = 0.019, fh = 0.197 and f to t — 0.283. The case of HE can be 
trivially added. Summation over % takes place in the case of mass-degenerate neutrinos. For 
the target neutrino density we shall use one helicity density with zero chemical potential, 
n Ui = 56 cm" 3 , corresponding to the temperature T = 2.73K of the microwave radiation, 
and £ = cos 2 9, sin 2 6 or 1, as explained in the footnote. 



Production of Z occurs through the interaction of flavor neutrinos, e.g. v e L + v e L in the case 
of Dirac neutrinos, or v e L + v^ L in the case of Majorana neutrinos. In practice one considers 
interaction of HE flavor neutrino, e.g., P e , with a physical mass-eigenstate target neutrino i>\ of 
mass mi. The probability to find this neutrino as v e is equal to cos 2 9 (or sin 2 9), where 9 is the 
mixing angle. Therefore, the cross-section in Eq.(|3lD must include £ = cos 2 9 (or sin 2 #). In the 
case of mass-degenerate neutrinos, m\ ~ m<i = m u , an incident HE v e interacts with v e component 
of f2 almost at the same resonant energy and £ = 1. There is no difference in the rates for Dirac 
and Majorana neutrinos: this becomes obvious if in counting the number of neutrino species one 
includes both vi and z^£ in the case of Majorana neutrinos. 
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Note that Eqs. (|32|) .(|33|) are exact formulae because integration in Eq.(|32|) takes place 
over a narrow resonant peak. 

The rate of Z° production with subsequent decay of Z° to an arbitrary channel / is 
given by riz(Z — > /) = nzff/ftot- Taking into account only dominant hadron channels, with 
fh/ftot = 0.696, and using the fact that pions transfer to electro-magnetic cascade half of 
their energy, one can find the energy density of electromagnetic cascade as 

uJcas = T^^-Eohz^to = 2irCa t n u .t I Vi {E )E^-^. (34) 

^ J tot J tot 

Using u cas < u° bs we obtain an upper limit 



ircr t n Vi t <m 4 z £, f h 



UEo) < V^47^ ( 35 ) 



One can see from Eq.(|3"q) that unless the production spectrum of neutrinos has a cutoff 
at some energy lower than E , the upper limit is provided by the lightest relic neutrino z/j. 
This is not surprising: in QCD spectra most of the energy is concentrated in particles of the 
highest energies. The lower m v is, the higher is E and the more energy it tranferred to the 
e-m cascade. 

The scale of neutrino masses suggested by oscillation solutions to the atmospheric 



neutrino and solar neutrino problems is m ~ V Am 2 . This gives the following masses: 
m u ph 5 ■ 10~ 2 eV from atmospheric neutrino oscillations, m u ~ 2 ■ 10~ 3 eV from SMA M SW, 
m u pa 4 • 10~ 3 eV from LMA MSW and m v pa 9 • 10~ 6 eV from VO solution. The latter mass 
requires a value of mx which is too large, so we disregard this case. 

It is easy to verify that we did not exceed the e-m cascade limit in the calculations of 
section IV. Indeed, Eqs.(|26|) and ([34]) with r 2 fj,/m x = 0.1, mx = 1 ■ 10 16 GeV, n v . = 
56 cm -3 , with Eq given by Eq.([3~0"l) and with both channels vi + z/j and vi + Vi taken into 
account, result in u cas = 1.2 ■ 10~ 7 / ^Jm\ eV/cm 3 , where mi = m u /l eV. The energy density 
u cas is well below the allowed limit for mi = 1 and is marginally below it for mi = 3 • 10~ 3 . 
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VI. UHECR FROM RESONANT NEUTRINOS 



Here we shall estimate the Ultra High Energy Cosmic Ray (UHECR) fluxes produced 
by neutrinos from hidden-sector TD. If the target neutrino density is enhanced in nearby 
large-scale structures (LSS), such as the halo (h) of our Galaxy, the Local Group (LG) and 
the Local Supercluster (LS), the large flux of the observed UHECR could be generated there. 
Photon fluxes dominate in the production spectra. The ratio of photon fluxes from a large- 
scale structure, I^ SS (E), and from extragalactic space, I* xtr (E), can be expressed in terms 
of the overdensity of target neutrinos in the structure, 5lss, and the length of gamma-ray 
absorption in extragalactic space, Ry(E), as 

I^ ss (E)/lf r (E) = 6 lss Rlss/R,(E), (36) 

where Rlss is the linear size of the large-scale structure. 

The overdensity factors for the galactic halo (h), Local Group (LG), and Local Super- 
cluster (LS) are discussed in the accompanying paper and for non-degenerate neutrinos 
they are estimated as 

b h v < 37ml b L v G < 13ml 5^ s ~ 1. (37) 

Here, ra\ = m„/l eV is the neutrino mass in units of 1 eV and the overdensity is defined 
as the ratio of the neutrino density with flavor i (antineutrinos are not included) in the 
structure to the average density of the same neutrinos in the Universe, n Vi = 56 cm~ 3 . From 
Eqs.(|36|),(|37|) one can see that while LS does not give an enhancement of UHECR flux, both 
the galactic halo (Rh ~ 100 kpc), and the Local Group (Rlg ~ 1 Mpc) give an enhancement 
of order (0.3 — l)mf. Note that this excess flux arrives without absorption. Estimates for 
both structures are given as upper limits, with the limit for the galactic halo being more 
reliable. In the estimates below, the index LSS refers to one of these two structures. 

Once again, let us take HS necklaces as an example of neutrino sources, with the diffuse 
neutrino flux I Vi (E) given by Eq . (|27|) . Using the formalism developed in Section V, one can 
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write down the flux of Z°-bosons with resonant energy E : 



I z o = 2£a t 8 1/ n Vi R LSS I n (E )E . 



Assuming a power-law spectrum of hadrons in Z° — > hadrons decay and using Eq. (|27|) for 
I Vi (Eo), one can easily calculate the UHE photon flux from the halo or LG: 




where Rlss is ~ 1 Mpc is the case of LG, and ~ 100 kpc in the case of galactic halo, 



where T had /T tot is the ratio of Z° decay widths, equal to 0.7. For p = 1.5 fc 7 = 0.116 and 
for p = 1.3 k 7 = 0.58. 

With 5 U = 5™ ax ~ ml, the UHE gamma-ray flux given by Eq.(^) is proportional to 
w? v and m x . For a fixed overdensity, the flux does not depend on the neutrino mass and 
depends only on the mass of X-particle as m x . 

As a numerical example let us consider the case of a gamma-ray flux from LG with 
two neutrino flavors and with a degenerate mass m u = 2 eV (£ = 1), taking the mass of 
X-particle mx = 1 ■ 10 15 GeV and £ = ( max = 10. For p = 1.5 and p = 1.3, the values 
of E 3 I^(E) at E = 1 ■ 10 20 eV are equal to 2.3 ■ 10 24 eV 2 /m 2 ssr and 1.8 ■ 10 24 eV 2 /m 2 ssr, 
respectively, i.e. close to the observed values. 

It is interesting to derive an upper limit for the UHE gamma-ray flux inside a LSS, using 
e-m cascade production in the space outside it. For LSS with a linear size Rlss an d neutrino 
overdensity 5 U , one obtains 



From Eq . (flUD one can see that the upper limit does not depend on mj and is proportional 
to ml~ p for 5 V = 5™ ax . For the parameters of LG (Rlg — 1 Mpc and S^ G = 13mf) 
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C = r 2 /j,/m 2 x , £ = cos 




(39) 




(40) 



and u cas = 1 ■ 10~ 6 eV/cm 3 one obtains, at E — 1 • 10 20 eV, fluxes equal to P 3 / 7 (P) = 
5.0 ■ 10 23 mf 5 eV 2 /m 2 ssr and 2.4 • 10 23 m5' 7 e^ 2 /m 2 ssr for p = 1.5 and p = 1.3, respectively. 
At m v > 2 both upper limits are consistent with observations. 

Turning the argument around, one can obtain a lower limit on the neutrino mass from 
the condition I™ ax (E) > I obs {E) at E = 1 • 10 20 eV: m v > 2 eV. 

More accurate calculations with realistic QCD spectra from Z° decay are given in the 
accompanying paper |HJ 

VII. DISCUSSION AND CONCLUSIONS 

Mirror matter is a natural option in models with G x G' symmetry, in particular in 
superstring models Eg x E' 8 . The coupling constants, the Higgs vev's, and the symmetry 
breaking patterns in the two sectors may or may not be the same, depending on whether or 
not the discrete P-symmetry interchanging the sectors is spontaneously broken. 

We assume that the two sectors communicate due to a non-renormalizable interaction (|l|), 
where the case A ~ m p i corresponds to gravitational interaction. These interactions result 
in neutrino masses and neutrino oscillations, including the oscillations between ordinary and 
mirror (sterile) neutrinos. 

Cosmological restrictions rule out a wide class of hidden-sector models; they are par- 
ticularly severe for the symmetric models in which the P-symmetry is unbroken. In such 
models, the number of light particles is doubled, and this is excluded by the cosmological 
nucleosynthesis. The nucleosynthesis constraints can be avoided by suppressing the temper- 
ature of the mirror matter, accompanied by an introduction of a lepton asymmetry, which 
suppresses v — > i/' oscillation. 

Inflationary scenarios resulting in different temperatures in the ordinary and mirror sec- 
tors can easily be constructed for asymmetric HS models. In the case of symmetric models, 
we discussed a two-inflaton scenario, first outlined in 0]. The inflatons and 0' belong to the 
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ordinary and hidden sector, respectively. In regions of space where 0' reaches the minimum 
of its potential earlier than 0, the products of 0' decay are diluted by the expansion driven by 
the ordinary inflaton 0. When also rolls to the bottom of its potential, we get a superhori- 
zon region dominated by ordinary matter. In stochastic inflation, superbubbles dominated 
by mirror matter are equaly often produced. We have shown that this scenario can work 
only with a suitable choice of the inflaton potential: for some choices the slow rolls of the 
two inflatons get synchronized, resulting in equal densities of mirror and ordinary matter. 
We also suggested an alternative version of the two-inflaton model where the potential has a 
metastable minimum. Then the inflating false vacuum decays through nucleation of 0- and 
0'-bubbles, and the segregation of ordinary and mirror matter is achieved in a natural way. 

Despite the suppression of mirror matter, hidden-sector topological defects can dominate 
over the ordinary ones. Once again, this can be easily arranged in asymmetric HS models. In 
symmetric models, the two sectors have the same types of defects with identical properties, 
but the cosmological densities of the defects need not be the same. We illustrated this pos- 
sibility by a two-inflaton model with a curvature-driven phase transition (see section III). In 
this model, HS topological defects are produced in a phase transition during inflation, when 
the mirror inflaton 0' is already at the minimum of its potential. The phase transition is 
triggered when the spacetime curvature (which is driven by the ordinary inflaton potential) 
decreases to some critical value. If this happens sufficiently close to the end of inflation, the 
resulting defects are not inflated away. The corresponding phase transition in the ordinary 
matter occurs much earlier, and ordinary topological defects are completely diluted by infla- 
tion. Thus, in the two-inflaton scenario we can have a desirable situation when the universe 
is dominated by ordinary matter and hidden-sector topological defects. 

HS topological defects produce high-energy neutrinos in the chain of decays of superheavy 
particles - constituent fields of the defects. All decay products are invisible in the ordinary 
matter and only mirror neutrinos oscillate into the ordinary world. The flux of neutrinos 
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from ordinary topological defects is limited by cascade photons which are produced in the 
same decays of pions as neutrinos. This restriction is absent in the case of HS defects. 
However, the cascade limit for mirror neutrinos reappears, though in a weaker form. After 
v' — > v oscillation, ordinary neutrinos produce hadrons, e + e~ and t + t~ in the 

resonant scattering off the background (dark matter) neutrinos: v + Vdm ~ > Z° ~^ hadrons, 
or . These particles (or products of their decay) initiate electromagnetic cascades on the 
microwave photons. The smaller is the mass of DM neutrinos, the stronger is the cascade 
upper limit [see Eq.([J^)]. This is because the resonant energy is inversely proportional to 
the neutrino mass: Eq = m 2 z /2m v . In QCD spectra, most of the energy is carried by high 
energy particles, and thus more energy is transferred to the e-m cascade when Eq is large. 

As a specific example of mirror topological defects, sources of high energy neutrinos, 
we studied the necklaces - magnetic monopoles connected by strings, with each monopole 
being attached to two strings. We found that, for a reasonable choice of model parameters, 
the diffuse neutrino flux can be three orders of magnitude higher than that from ordinary 
necklaces, being still consistent with the cascade upper limit imposed by the resonant pro- 
duction of Z° bosons. Note, however, that the accuracy of our calculations is limited by the 
power-law approximation of the energy spectra. 

A diffuse flux of UHE neutrinos from HS topological defects can produce the observed 
flux of UHE cosmic rays due to resonant interaction with the Dark Matter neutrinos in the 
Local Group, if the mass of the target neutrino is m u > 2 eV. The atmospheric-shower 
producing particles in this case are UHE photons. Their spectrum does not exibit a cutoff, 
because of the relatively small size (R ~ 1 Mpc) of LG. 

This model of UHE cosmic rays requires mass degenerate neutrinos with m u > 2 eV. Our 
model for neutrino masses and mixing is described in the Appendix. Mirror muon neutrino 
v' is maximally mixed with the ordinary neutrino, and both have masses m v ~ 2 eV. 
Their mass difference, Am 2 pa 2 • 10~ 3 eV 2 , is responsible for the atmospheric neutrino 
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oscillations observed in Super-Kamiokande. Solar neutrino anomaly is explained by v e — > 
LMA MSW solution with Am 2 w 4 ■ 10~ 5 e\/ 2 and sin 2 29 « 0.80. Thus, we assume that the 
three neutrinos, i/ e , z/ M and i/' are maximally mixed and mass degenerate (m v ~ 2 eV). 

The calculated neutrino flux is below the upper limits obtained from horizontal air shower 
observations at EAS TOP @ and AG ASA Q at 10 6 -10 7 GeV and marginally below Fly's 



Eye limit |35] at 10 11 GeV. The predicted neutrino fluxes can be detected by this technique 



with the help of these and future bigger arrays, like e.g. "Auger" detector |36]|. However, 



the best hope for detecting these neutrinos probably rests with the future satellite detectors 



such as OWL (Orbiting Wide Field Light Collector) and AIRWATCH |8 
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APPENDIX A: MASS DEGENERATE NEUTRINOS IN MIRROR MODELS 



We consider G x G' model with EW symmetry SU 2 (L) xUiinG and SU' 2 (R) x U[ in G' . 
G and G' representations communicate through operators of dimension d = 5 with a scale 
A < m pl . 

The particle content of the EW group relevant to the neutrino masses is 



R, 



V 



-0;/ 



(Al) 



for SU 2 (L) x Ui, and 
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(A2) 



for SU' 2 (R) x U[. Here, and 0' are the Higgs fields and I and V are charged leptons. 

There are no light singlets Vr and v' L in our model. There are 4 neutrino states: 
ul, z/£, i/jj, v'fi, in terms of which the most general expression for the mass matrix is 



C ~ ( 



m L M 
M m R 



\ 



(.A 



+ h.c. 



(A3) 



Since u' ' R v c L 



vlv'r, there are three independent mass operators in Eq. (|A"3"|) : 



vlVri v l ij li v> r v 'ri an d they are generated in our model. Indeed, these operators are: 



SU 2 (L) x Z7i : (W c ) (0 C ^ 
Sl£(J2) x U[ : (^0 /c ) (0 /c ^ 
intergroup : (V>l0 c ) \ d>' cJr ^ 



vlVr 



(A4) 
(A5) 
(A6) 



Arrows show the neutrino mass operators after the EW symmetry breaking, < 0o >= vew 
and < <p' >= v' EW . Communication of visible and mirror sectors is accomplished by the 
intergroup term in Eq . (|A6|) . It has a dimensional scale A, and thus one obtains 



M = vewv' ew /A. 



(A7) 



This is the basic neutrino mass scale in our model, and we want it to be M ~ 1 eV. For 
v' ew /vew ~ 10, we need A ~ 10 14 GeV. 

The scales of d=5 terms operating inside SU2(L) xU\ and SU' 2 (R) x U[ groups, A^ and 
A R , can be different. For our model (see below) we need the following hierachy of masses: 



rriL < rriR <C M. 



(A8) 



It can be provided by A Ar < A E , or in the model-dependent way. One can observe, for 
example, that both intragroup terms ( |A4j ) and (|A5|) violate the lepton number defined for 
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the doublets as = L$> = 1, while the intergroup operator (|A6|) conserves it. One can 
build a model with one universal A where the intragroup d=5 operators are forbidden and 
thus m L and m R are suppressed. 

Let us assume a local U\ symmetry for massless particles before symmetry breaking, with 
the following charge assignment: q = +1 for ipL, Ir, iP'r, I'l an d g = for 0, 0'. 

The terms ( |A4| ) and (X5) do not conserve q. Let us introduce two new scalar SU2 singlets, 



$ and with q = 2. Now, apart from the operator (|A6|) , we can write two other SU2- 
singlet operators conserving the charge q and the electric charge. They are the following 
d = 6 operators: 

(M c ) (0 c ^£) , ^ (^0' c ) («) . (A9) 

After EW and Ui symmetry breaking with vev's < $ >= V and < $' >= V', respectively, 
one obtains 

M = v EW v' EW /A, m R = v'£ w V'/A 2 } m L = v 2 EW V/A 2 } (A10) 



which satisfy the hierarchy ( A"8 



Let us now come back to the mass matrix (|A3| ) . Its diagonalization gives the masses of 
eigenstates and the mixing angle for visible and mirror neutrinos: 



mi, 2 = i^rriR + m L ± y 4M 2 + (m R - m L ) 2 j w 2M (All) 

Am 2 = TOj - mj ps 2(m_R + m L )M (A12) 

sin 20 = 2M ps 1 (A13) 
/4M 2 + (mjj + m L ) 2 



The hierarchy of masses ( |A^j|) provides mass degenerate neutrinos with (almost) maximal 
mixing. 

Till now we considered asymmetric mirror models. In the case of symmetric models 



12] , the masses M, mi and m R are considered as free parameters, and thus the hierarchy 



condition (|A8|) can be arbitrarily fulfilled. 
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As a realistic example we can consider the case when — > u' oscillations explain the 
atmospheric neutrino anomaly (Am 2 « 2 ■ 10~ 3 eV 2 ) and LMA MSW explains the solar 
neutrino problem [y e — > oscillation). In this case, 2M ~ 1 eV and m^ ~ 1 • 10~ 3 eV. 
All three neutrinos, z/ M , v e and z/^, must have mass 2M ~ 1 eV and be bimaximally mixed. 
The oscillation length (z/ — > z/ M ) for HE neutrinos with resonant energy Eq is only l osc = 
AnE / Am 2 = 5-10 20 cm. The values of Am 2 in this case are 2-10 -3 eV 2 and ~ 4-10 -5 eV 2 for 
the atmospheric and solar neutrinos, respectively; they are not affected much by radiative 
corrections [j^ . Radiative splitting affects mostly r neutrino, because of the large Higgs 
coupling with the tau. This effect is not important in our scenario where only three other 
neutrinos, z/ M , v' and v T , are degenerate in mass. 
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